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Abstract 

There has been much recent work on the spectrum of the random 
non-hermitean Hamihonian which models the physics of vortex line pinning in 
superconductors. This note is loosely based on the talk I gave at the conference 
"New Directions in Statistical Physics" held in Taipei, August 1997. We describe 
here new results in spatial dimensions higher than one. We also give an expression 
for the spectrum within the WKB approximation. 



*To appear in the proceedings of "New Directions in Statistical Physics," Taipei, August 
1997, Physica A, to be published. 



1 Introduction 



A non-hermitean Hamiltonian inspired by the problem of vortex line pinning in su- 
perconductors has attracted the attention of a number of authors ^ ^, ^, ^, 
|, |, 0. In one version, we are to study the eigenvalue problem 

Hijijj = ^[e^iJi+i + e-^ i/ji^i) + WiiJi = Ei)i, i,j = 1,...,N (1) 

j 

with the periodic identification i + N = i of site indices. Here the real numbers zfj's, 
the site energies, are drawn independently from some probability distribution P (w) 
(which we will henceforth take to be even for simplicity.) This Hamiltonian describes 
a particle hopping on a ring, with its clockwise hopping amplitude different from its 
counter-clockwise hopping amplitude. On each site there is a random potential. The 
number of sites is understood to be tending to infinity. We will also take t and h to 
be positive for definiteness. (The hopping amplitude t can be scaled to 2 for instance 
but we will keep it for later convenience.) 

Note that the Hamiltonian is non-hermitean for h non-zero and thus has complex 
eigenvalues. It is represented by a real non-symmetric matrix, with the reality im- 
plying that if E is an eigenvalue, then E* is also an eigenvalue. For P (w) even, for 
each particular realization {wi} of the random site energies, the realization {—Wi} 
is equally likely to occur, and thus the spectrum of H averaged with P (w) has the 
additional symmetry E — E. 

Without non-hermiticity (/i = 0), all eigenvalues are of course real, and Anderson 
and collaborators |T^ showed that all states are localized. Without impurities {wi = 
0), Bloch told us that the Hamiltonian is immediately solvable with the eigenvalues 

E{6) = t cos {6 — ih) = cos 9 cosh h + i sin 6 sinh h , (2) 

with 6 = (n = 0, 1, ■ ■ ■ , A^ — 1), tracing out an ellipse. The corresponding wave 

functions ^j""* ~ e*^-' are obviously extended. We are to study what happens in the 
presence of both non-hermiticity and the impurities. 
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Figure 1: The eigenvalues of one particular realization of the Hamiltonian in (1) with 
= 500, t = 1, h = 0.4 and with the random site energies Wi equal to ±0.5. 



In figure (1) we show a typical spectrum for the Hamiltonian in (|l|). The spectrum 
is made of two parts, a more or less elliptical (or "lemon-like") curve, consisting of 
complex eigenvalues, and "two wings," a right wing and a left wing, consisting of real 
eigenvalues, poking out of the complex part of the spectrum. 

In the literature [|l], |^, [Tl|, [1^, it was argued that localized states have real 
eigenvalues and that states of H corresponding to complex eigenvalues are extended, 
that is, delocalized. Remarkably, non-hermitean localization theory is simpler in this 
respect than the standard hermitean localization theory of Anderson and others |14 |. 
To understand the localization transition, we have to study only the energy spectrum 
of H. 

1 refer the reader to the talk by N. Hatano in these proceedings for an illuminating 
discussion of this subject. 

In the talk at the conference, I described work done with J. Feinberg 0, |l^ and 
with E. Brezin [||. Since the material covered can be found discussed in detail in these 
references, it may be more useful to the reader for me to describe here some slight and 
unpublished work done with E. Brezin. We will mention two separate topics here. 
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2 Continuum Limit and WKB Approximation 



As is obvious and as has been discussed in the hterature, we can easily take the 
continuum hmit of (|I|) by letting x = ia, i oo, a 0, t oo, to? = — ^, 
e = E — t, and A = ^. Note that we have introduced the gauge potential A, equal 
to h divided by the lattice spacing, and hence with the usual dimension of inverse 
length. We obtain 

- 7^ f(V + Af - A^] ij + Vifj = ei/j (3) 
Inserting ijj = e*'^ we have (setting 2m = 1) 

S'^ - 2iAS' - iS" = e - V (4) 
The WKB approximation involves neglecting S", as usual. We obtain 



S' = lA ± ^J-A^ + {t - V) (5) 
with the boundary condition 

S{L) = S{0) + 2Tin (6) 

Let us define 



* - ^ (7) 
In the limit L oo, k becomes a continuous variable. Then we find 

(k - lA) = ± ^ dx Ve - A^ - V (8) 
L Jo 

One way to look at this expression is to study the effects of the non-hermitean 
potential A. In the absence of the non-hermit icity, we have 



We thus obtain a remarkably compact expression for the WKB spectrum of 



e{k) = eo{k - iA) + 



(10) 



Note that we have to know the usual hermitean spectrum £0(9) for complex values of 
its argument. In particular, we can study the effect of non-hermiticity perturbatively 



into the complex plane. For instance, if eo{k) = eo + ak"^ + bk* + . . ., then 
e(k) = Co + ^2(1- a) - 2ia Ak + ak^ + .... 

This discussion holds whether V is random or not. For random V, it would 
be interesting to see if it would be possible to obtain a relation between the non- 
hermitean density of states 

p{x,y) = { J dkS(x - A^ - Re eoik - iA)) 6 {y - Im to {k - iA))) (12) 
and the hermitean density of states 





(13) 
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3 Higher Dimensions 

Much of the hterature on this subject is numerical and, with some exceptions, focused 
on one dimensional space. In real systems, vortices move in two dimensional space. In 
a recent work [|15| , Nelson and Shnerb studied the higher dimensional problem using 



numerical, perturbative, renormalization group, and other techniques. Thus, it may 
be of some interest to report here some exact analytic results in higher dimensions. 

First, we have to review some basic formalism. As is standard, we are to study 
the Green's function 



G{z) = (14) 
z — H 

The averaged density of eigenvalues, defined by 

p{x, y) = {j^T. <^(^ - S{y - ImE,)) , (15) 

i 

is then obtained (upon recalling the identity ^ ^ = tt6{x) S{y)) as 

1 d 

<{x,y) = 

with z = X + iy. Here we define 



TT OZ 



G (z) ^ {- tr G (z)) (17) 

and < . . . > denotes, as usual, averaging with respect to the probability ensemble 
from which H is drawn. It is also useful to define the "bare" Green's function 



G.(z) = ^ 



in the absence of impurities and the corresponding Go{z) = ^ ti Go{z). For the 
spectrum in @), we can determine Go{z) explicitly [§]. We need hardly remark that 



in determining Go{z) no averaging is involved. 
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In a formula relating the Green's function G to Gq was obtained. Here we give 
an alternative derivation, perhaps somewhat more direct, of this relation. From 



H = + Y.^kPk (19) 

k 

where we defined the projection operator 

Pk = \k){k\ (20) 
onto the site and the definitions of G and Gq we have 



= G^' - Y.WkPk = G^\l - GoWk Pk) (21) 

k k 

and hence 



G 



, 1 - Go Wk Pk I 
\ k / 



Go (22) 



We might consider expanding this expression as a power series in the WkS and then 
average with P{w). Indeed, that would be the correct procedure for many distribu- 
tions. However, for some P{w) such as the Cauchy distribution 

^H = -^^' (23) 

with its long tails extending to infinity, the moments {w^), for n an integer, do not 
exist. The physics here is that we have to first sum the effect of repeated scattering 
on the impurity potential at one particular site k 

1 - Wk{k \Go\ k) 1 - Wk Go 
The last equality follows from the assumed translation invariance of Ho- Since the 
probability distribution P{w) is of course normalizable, the moments (v^), for n an 
integer, in contrast to {wji), always exist. Note that in contrast to Wk, Vk is a function 
of z. We can solve for Wk in terms of Vk as 
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Vk 



Inserting (|25|) into 



G 



1 + Vk Go 
we obtain without further ado 



(25) 



/ VkGoPk 
V \l + VkGoPk 



n -1 



Gq 



VkGoPk 

^l + VkG'o 



1 -1 



Go (26) 



The reader is invited to average the trace of ( PB| ) with his or her favorite P (w), and 
thence to obtain the density of eigenvalues p{x,y) by (|1^). 

In general, it is non-trivial to carry out this averaging explicitly. For the Cauchy 
distribution, however, we have the drastic simplification that 



with 



I 7 



1 - Z7 Go ' 



for Im 2; < 



and 



(27) 



(2^ 



-I 7 



= , for Im 2; > 



1 + 27 Go 

which by translation invariance is independent of k. Defining as usual 



(29) 



we see that for the Cauchy distribution 



(30) 



{v) 



1 + {v} Gq 

and is thus equal to =FZ7 in the upper and lower half plane respectively. We thus 
obtain from ([TT 



G{z) = Go(2-S) 



(32) 
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Contrary to some erroneous claim in the literature, for an arbitrary probability 
distribution T,{z) is in general a matrix and (^) does not hold. 

For the Cauchy distribution, however, we obtain from (PT|) and (ET) 



G{z) = Go{z + 27) e (Im z) + Go{z- 27) 9 (-Im z) (33) 

The density of eigenvalues are given by ([16|) . When ^ in ([16|) acts on Go in (|33|) 
we obtain the "central part" of the spectrum 

Pcentraiix.y) = po{x, y + 9{y) + po(a;,l/ - 7) ^(-y) (34) 

In other words, the random site energies push the density of eigenvalues po{x,y) of 
the non-random problem towards the real axis by a distance 7. This accounts for 
the complex eigenvalues. In figure (1) this central part of the spectrum consists of 
a curve |^, |10]. However, in general, in higher dimensional space for instance, the 
central part of the spectrum forms a two dimensional blob. 

When ^ in (|I^) acts on the step functions in (|53| ) we obtain the density on the 
two "wings": 



7 1 1 
tr 



Pwing{x,y) = -7^^5(2/) (Go{x + i-f) - Go{x-i-f)) = 6{y) u \2 , 2 

ZTTl ^ ^ TT iV [X — SIq)^ + 7^ 

(35) 

We thus have the general result that Pmng is constructed out of the sum of Cauchy 
distributions centered about each of the eigenvalues of Hq. At first sight, this result 
seems puzzling, since Pwing is not obviously real. The resolution is that we have to 
look more closely at how the averaging over the Cauchy distribution is done. The 



discussion in the literature [|1^] is given for hermitean Hq^s, and hence had to be 
generalized to non-hermitean Hq^s in P). Imagine expanding ( p6D in a series in the 
VkS. According to (^71) we have to average v = -^jz^j^, that is, to integrate vP{w) 
over w. In [^, Gq{z) is explicitly known, and we can see that for z in the upper 
half plane, the integration contours in w can be closed in the lower half plane and by 
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Cauchy's theorem w is effectively set equal to — ^7, and vice versa if z is in the lower 
half plane. 

In the present discussion, we have to find a class of Hq^s for which this procedure 
continues to hold. It is not too difficult to find such a class. Suppose the eigenvalue 
spectrum of Hq is such that if E is an eigenvalue, then E* is also an eigenvalue. As 
was mentioned in the introduction, this is the case if Hq is real (but not necessarily 
symmetric.) Then Go{z) = Go{z*) and it is straightforward to show that the 
integration procedure outlined above holds. It is an interesting question whether 
there exists a class of Hq^s broader than the one found here for which this integration 
procedure still holds. 

We now see, gratifyingly enough, that for the class of i^o's found here Pwing is 
indeed real and positive definite. Indeed, we can write (|35|) entirely in terms of 



Pwing{x,y) = 6{y) - / dudvpo{u,v)- 212 , ^ 2^ \^ (^6) 

where we used po{u,v) = po{u, —v). 

Thus, we simply have to evaluate the Green's function for a given Hq, plug in 
(|3^) and (pSf), and obtain the density of states for H. We emphasize that in this 
derivation the only thing we assumed for Hq is its translation invariance. Our result 
applies even if Hq contains non-nearest neighbor hopping, for instance. In particular, 
this result holds for any spatial dimension. 

As an example, in we applied this result to the Hq in (|I|) and obtained 



1 2 2 Ji^-x^ + f^ + B{x,-f,t) 
Pwing{.x,y) = -j=- 6{y) e{x -x„i„(7)) ^7 7\ (37) 

V ZTT ID[X,'J,l) 

where we have defined 

B{x,j,t) = ^(^2 + 72)2 + 2t2(72_2.2) + f4 (38) 

and 
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The critical value of 7 at which the wings extend over the entire real axis is determined 

by Xminil) = 0: 

jc = t sinh h (40) 

For randomness larger than 7c, all states are localized. 

Note that the density of eigenvalues on the wings does not depend on the value of 
the non-hermiticity h at all; h comes in only in determining where the wings end. This 
explicitly verifies the argument [|1], 0, |Tl], 0, |13[ in the literature that the localized 
states do not know about the non-hermiticity. 

We now understand that most of the qualitative features exhibited by these results 
obtained in p are generic for the Cauchy distribution. From (^), we see that 7ci, 
defined as the critical randomness at which wings first appear, is equal to O''" for 
any translation invariant Hq. We also obtain a general result for 7c, the critical 
randomness at which all states are localized, namely 

7c = Umax (41) 

where ymax is the maximum value of y in the support of po{x,y). 

In summary, in any dimensions, for any Hq whose eigenvalues come in complex 
conjugate pairs, and with the random site energies governed by the Cauchy distribu- 
tion, we have the result that the density of eigenvalues is given by the sum of two 
terms, Pcentrai{x,y) iu (|3|), and Pwing{x,y) in (^ which has the form 6{y) f{x). By 
the way, this explicit delta funtion accounts for the "ridge" along the real axis seen in 
the numerical work of Hatano and Nelson (see figure (19c) in Q]). In the numerical 
result, we see a background consisting of a two dimensional blob on which eigenvalues 
are more or less uniformly distributed. Superimposed on this background is a dense 
cluster of eigenvalues on the real axis, the so called "ridge." (The entire spectrum is 
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thus reminiscent of a cross section of a galaxy, with a halo and a concentration of 
stars in the galactic plane.) 

The numerical work of Hatano and Nelson was actually performed for the box 
distribution P{w) = ^diV"^ — w"^). We thus have an indication that some features 
of the spectrum are generic. One feature which can clearly be attributed to the fact 
that the Cauchy distribution has a long tail is that our Pwing{x,y) in (^) extends to 
X = ±00, falling off as Invoking the galactic analogy, we can describe the situation 
for Cauchy distribution by saying that the halo is concentrated near the center of the 
galactic plane, with the density of stars thinning out rapidly, while for distribution 
with finite support the halo more or less fits over the galactic plane. 

A quantity of some interest is the fraction of eigenvalues in the "galactic plane", 
which we can obtain by integrating Pwingi^, y). We find 

/oo /■ 
pwingix,y) = 6{y) / dudvpo{u,v)9{--f < i; < 7) (42) 
-00 J 

00 

Note that / dx pwmg {x, y) = S{y) for 7 larger than 7c, in accordance with (p]). 

—00 

The reader is invited to evaluate Gq{z) for his or her favorite Hq. 

We will give two illustrative examples here. As a first simple example, we can 
consider a ladder, on whose two legs are two copies of the Hamiltonians in (^, and 
with hopping between the two legs given by the amplitude t'. More precisely, we take 



Hq ai, f3j = I (e'* 6ai3 (^i+l, j + ^ ^ ^a/B 

+ I Sij {6a+l, 13 + Sa, l3+l) (43) 

where a, (3 take on two values and denote the two legs. The periodic boundary 
condition i = i + N is imposed as before. Alternatively, we can think of a particle 
carrying a binary label (spin, flavor, etc) hopping on a ring. This model was studied 
numerically by Hatano. See these proceedings. 
The eigenvalues of Hq are obtained immediately 

f 

E{e) = t cos {9 - ih) ± - (44) 
11 



tracing out two overlapping ellipses separated by t' . The corresponding Go is the sum 
of two pieces 



Go{z) = l^z + + go - (45) 

where gQ{z) is the Green's function for the spectrum in @) 

Now we turn on Cauchy randomness for the site energies. According to (p3) and 



(|35|) , the spectrum is just the superposition two "squashed ellipses," each with its 
own wings poking out. For comparison, we show the numerical result in figure (2). 

We will finally venture into two dimensional space. We can of course immediately 
write down the analog of (0) for a lattice of any dimension. For two dimensions, we 
denote the parameters for hopping in the x direction by t, h and for hopping in the y 
direction by t\h' . The spectrum of Hq is obviously given by 

E{e, 9') = t COS {9 - ih) + t' cos {9' - ih') , (46) 

with 9 = 9' = (n, n' = 0, 1, ■ ■ ■ , A^— 1). The spectrum is a two dimensional 

blob, constructed geometrically by placing an ellipse centered on each and every point 
of another ellipse. See figure (3). 

Depending on the values of t, t', h, and h', the spectrum can have different topolo- 
gies. Let us refer to the case in which t = t' as the isotopic hopping case. Note 
that in that case the "hole" in the spectrum closes when the hermiticity field {h, h') 
points in the (±1,±1) direction. (It is straightforward to take the continuum limit. 
In general, the mass m in (|^) becomes anisotropic; in the isotropic hopping case, we 
obtain (0) with the replacement A ^ A = ^(/i, /i').) 

The corresponding Green's function 

^ / ^ f d9 d9' 1 

can be obtained by quadrature. We would not bother to carry out this computation 
here. Instead, for illustration, we will compute this integral only in a simplifying 
limit. 
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Figure 2: The spectrum for the ladder problem (with = 200): (a) For t' large 
(here t' = 2 with t = 1 and /i = 0.8), tiki two ellipses are separated. The spectrum 
is shown for 7 = 0.3. (b) For smaller t' (here t' = 1, with t = 1, h = 0.8), the 
two ellipses overlap. The spectrum is shown for 7 = 0.1 (c) Same as (b) but with 
7 = 0.4 




ReE 



Figure 3: The spectrum of Hq in two dimensional space (ellipse + ellipse) for (a) 
t = 0.5, h = 0.3, t' = l,h' = 0.6 and for (b) t = 0.5, h = 0.6, t' = 1, h' = 0.2. 
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In 0, it was pointed out the the problem (|1]) simphfies in the maximally non- 
hermitean or "one way" limit, in which the parameters in (|l|) are allowed to tend to 
the (maximally non-hermitean) limit /i — > oo and t ^ such that 

(48) 

The particle in (|l|) can only hop one way. The spectrum (0) becomes 

E{9) = e'\ (49) 

with 6 = ^ , (n = 0, 1, ■ ■ ■ , — 1) and the ellipse associated with (0) expands into 
the unit circle. The corresponding Green's function can be determined immediately 



d9 1 



2tt z- e^^ 
9i\z\'-l) 



(50) 



z 

(We trust the reader not to confuse the Heaviside step function with the angle vari- 
able.) For I z I < 1, Go{z) = 0, and for | z | > 1, Go{z) = \. By (0) the 
density of eigenvalues vanishes in both regimes. The density of eigenvalues indeed 
has support on the unit circle, in accordance with (^). 

In two dimensional space, we can obviously take either the single "one way" limit 
in which hopping along the y direction becomes one way, or the double "one way" 
limit in which hopping along both the x and y directions becomes one way. From 
(HBI), we see that the spectrum changes from "ellipse + ellipse" to "ellipse + circle" 
to "circle + circle." We show the spectrum for the single "one way" limit in figure 
(4) and for the double "one way" limit in figure(5). 

The double "one way" limit is taken by letting h,h' ^ oo and t, t' 0, but in 
general the limit can be taken such that ^ ^ ^ e^' . For the sake of simplicity, we 



2 ^ ' 2 

t h _ i_ „h' 
2 ^ 2 



will take ^ = ^ e'^ ^ 1. It will be obvious to the reader that the Go(z) for the 



more general anisotropic case can be obtained from the expressions below by a trivial 
rescaling. 
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Figure 4: The spectrum of Hq in the single "one way" hmit for t = 0.5, h = 0.3 
(elhpse plus circle). 
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Figure 5: The spectrum of Hq in the double "one way" hmit (circle plus circle). 
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So, "just for fun," we will evaluate 



de_ dff_ 



27r 27r z - e^^ - e^^' 



z — e 



ie 12 



2n 



(51) 



Similarly to our preceding discussion for the one dimensional case, for \ z \ > 2, 
Gq{z) = ^, and by ([T6|) the density of eigenvalues vanishes there. By rotational 
invariance, we can take z to be real and equal to x such that < a; < 2, and obtain 



Goix) 



1 



2 TT X 



de 



cos 9 — X 



cos 9 



2x 



After some straightforward computation, we find for | 2; | < 2 



1 r 
— arc cos - 

TT 2 



— arc tan 

TT 



1 — r 

1 + r 



2-r 



where r = \ z \. The density of eigenvalues is given by (|T6D to be 



(52) 



(53) 



Poir) 



1 dGn 



TT dz* 

1 1 du 

4iT r dr 
1 1 



(54) 



where u{r) is the function defined by the square bracket in (|53D. Note that poir) 
diverges as r — 2 from below. The density of eigenvalues can of course also be 
obtained directly from its definition 

{x,y) = { j d9S{x - Re {e'' + e''')) 6 [y - Im {e'' + e'''))) (55) 



When we turn on the Cauchy random site energies, the density of eigenvalues is 
then given by (|3^) and (|3^). Since we have po explicitly, we can work out various 
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quantities. For example, the fraction of eigenvalues on the real axis is, according to 
(p|), given by (for 7 < 2) 



4 ri rVA—v'^ '[ 

frealil) = ^ dv du (56) 

We can obtain the u integral in terms of an elliptic function and then determine 
frealil) numerically. We show in figure (6) a plot of freaiil)- 




Figure 6: The fraction of eigenvalues for freai as a function of randomness 7. 

We show some numerical results in figure(7) for two dimensional hopping with 
Cauchy randomness. Define the quantities youter = t' sinh h' + t sinh h and 
Uinner = t' siuh h' — t siuh h. FoT the parameters used in figure(7), we have 
y outer = 0.789 and y inner = 0.484. According to (H) and (|), youter and y inner 
determine the intersection of the outer and inner edges of the spectrum with the 
imaginary axis, for the range of the parameter values chosen here. (The complete 
determination of the outer and inner boundaries of the spectrum is straightforward 
but tedious, as it involves the solution of a quartic equation. We will content ourselves 
here with a remark about the simple case t = t' in which yinner = t{smh.h' — sinh/i). 
This expression is however not correct for all values of h; in particular, not for h = 0, 
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as we can see geometrically. For h small, we see from the geometrical construction 
"ellipse + ellipse" that the correct expression is Uinner = i(sinh h' tanh h'). The 
cross-over value of h is then determined by sinh h = sinh h' — sinh h' tanh h'.) 

We have drawn in the values of y outer and y inner for the values of t,t',h, and h' 
chosen here as the two horizontal lines in figure (7a); we see the outer and inner edges 
of the numerical spectrum indeed are equal to y outer and y inner, as it should be. As 



discussed before (JSJ), when we turn on the Cauchy randomness, the spectrum gets 
pushed towards the real axis by a distance 7. Thus, the critical value of 7 for the 
hole to disappear is given by 



7c, hole = t' sinh h' — t sinh h (57) 

(for values of t', h', t, and h such that a hole exists, and such that yinner is positive.) 
For figure (7), 7c, hole = 0.484. We see from figure (7b) that at 7 = 0.2 the pres- 
ence of the hole is clear. (The horizontal lines in figure (7b) represent {youter — 7) 
and {yinner — 7)0 In figure (7c), the hole has indeed disappeared. The fluctuations 
at finite can be studied using the density density correlation function defined in 0]. 
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Figure 7: The spectrum of if on a two-dimensional 20 x 20 lattice with t = 0.5, 
h = 0.3, t' = 1, h' = 0.6: (a) without randomness (compare with figure 3a), (b) 
with Cauchy randomness of strength 7 = 0.2, (c) with Cauchy randomness (b) with 
7 = 0.5. Note the resemblance to a galag|ic cross section, as mentioned in the text. 
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